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Abstract 

^| We construct intrinsic on-and off-diagonal upper and lower estimates for the transition prob- 

r— I ability density of a Levy process in small time. By intrinsic we mean that such estimates 

reflect the structure of the characteristic exponent of the process. The technique used in 
the paper relies on the asymptotic analysis of the inverse Fourier transform of the respective 
characteristic function. We provide several examples, in particular, with rather irregular Levy 
measure, to illustrate our results. 

t" — Keywords: transition probability density, transition density estimates, Levy process, Laplace 

! ! method. 

Oh 

g 1 Introduction 

A wide variety of asymptotic results on the small time behaviour of the transition probability 
^ density for Levy and, more generally, Levy driven and Levy-type Markov processes is available 
00 in the literature. In |Le87j and jlsh94j this problem was treated by means of the Malliavin 
calculus, see also |Pi97aj . [Pi97bJ, [IshOlJ, where the Picard's modification of the Malliavin calculus 
introduced in |Pi96] was used. Another approach, applicable for Levy-type Markov processes and 
developed in | BBCK09] . |CKK11|, [C09] , |BGK09j . |CK03] . |CK08| . |CKK09] . is based essentially 
O on the method introduced in [CKS87J, and involves the Dirichlet form technique. We refer also 
^ to more specific papers |StlOa] . |StlOb] . |RS10j . where the particular structure of tempered stable 
processes is used, and to |BNOO| . jFOllj . |FH09j . [RW02J for other results on the small time 
asymptotic expansions. Of course, this list of publications is far from being complete. 

Typically, the results available in the field involve some restrictions on the initial Levy noise, 



X 

C\3 formulated in an explicit and prescribed form, e.g., the two-sided power bounds for the truncated 



second moment of the Levy measure in |Pi97a] . |Pi97bj . and |Ish01| : see (2.10) below. On the 
contrary, the approach developed in the current paper is free from any assumption of such a type. 
We introduce two auxiliary functions ip L and ijj u , see (2.3) and (2.4) below, which give natural 



lower and upper bounds for the real part of the characteristic exponent tp of the Levy process, see 
(2.5) below. Our principal assumption is that the trivial inequality ip L < ip u is invertible, in a 



sense; see (2.6). This assumption is satisfied for a rather wide range of Levy processes including, 



for instance, the well studied stable and "a-stable like" processes, see Examples 4.1 and 4.2 On 
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the other hand, this assumption may hold true as well even when the real part of the characteristic 
exponent of the process exhibits a rather irregular asymptotic behavior at oo, see Example |4.3 



Under assumption (2.6 ) we are able to give two types of small time estimates for the distribution 
density pt(x) of the process Z t . First, we prove the on-diagonal estimates, i.e. two-sided bounds 
for the maximal value of the distribution density: for given to > there exist positive constants c 
and d such that 

cp t < maxp t (x) < dp t , t e (0,t Q ], (1.1) 



with 



p t =[neif>Y~ 1 \i/t) 



where 

/ \ hi] 

(Reip) (x) := inf{C > : Re^(C) > x}, x>0 

is the quasi-inverse function of the real part Keip of the characteristic exponent. Second, we give 
more informative, but more complicated, two-sided off-diagonal estimates. The structure of such 
estimates is a subject of separate discussion. It is well known (e.g. |Ish94] . |Pi97a| . [Pi97bJ, [IshOlJ), 
that for the transition probability density of a Levy driven Markov process such a structure, in 
general, should be more complicated than a "bell-like" estimate, which is typical for the diffusion 



case. We discuss this feature in details in Section 2.2 and in Examples |4.1[ [42] below. Motivated by 
this observation, we propose a new type of non-uniform estimates, which we call compound kernel 
estimates. For both types of the estimates described above we use the term intrinsic, because the 
structure of the estimates is mainly determined by the characteristic exponent and Levy measure 
of the process. In the main results any prescribed function, e.g. power-type, regularly varying, 
etc. is not involved neither in the assumptions, nor in the estimates themselves. In addition, we 
discuss some auxiliary assumptions which make it possible to simplify our main estimates, and to 
present them in a more compact form. 

A key ingredient in our approach is the asymptotic analysis of the inverse integral Fourier 
transform of the characteristic function of the initial process, which involves a proper version of 
the Laplace method (e.g. |Co65j ). and is a "small-time counterpart" of the asymptotic analysis 
developed in |KK11] and |K12] . For such an analysis, it is important to have the respective 
characteristic function explicitly. Consequently, the range of applications of such an approach 
seemingly does not exceed the class of Levy processes (or, slightly more general, Levy driven 
stochastic integrals with deterministic kernels). The next natural step in this research direction 
would be an extension of the asymptotic results obtained in this paper to Levy driven and Levy- 
type Markov processes using a perturbation technique, analogous to the parametrix method for 
the diffusion processes. This is a subject of a separate forthcoming paper |KK12] . 

To make the exposition as transparent as possible, we formulate our main estimates for one- 
dimensional Levy processes. Multidimensional generalization seemingly would not contain any 
substantially new difficulties. We postpone such a generalization to further publications. 

We now outline the rest of the paper. In Section [2] we introduce the notation and the assump- 
tions, and formulate our main results. Proofs are given in Section [3] In Section [4] we introduce 
and discuss some examples. 
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2 Main results 



2.1 Notation, assumptions, and on-diagonal bounds 

In the sequel, Z t is a real-valued Levy process with the characteristic exponent ip, i.e. 

Ee l ^ Zt = e~ m \ £ G R. 

We assume that Z t does not contain a Gaussian component. In this case if) admits the Levy- 
Khinchin representation 

V(0 = ia£ + [ (1 - ^ u + i£ul H<1 )n(du), (2.1) 

where a G R and /x is a Levy measure, i.e. / R (l A u 2 )fi(du) < oo. We will treat separately the 
situation where the process Z t is symmetric, which means in terms of the characteristic exponent 
that if) is real-valued. In what follows we assume 

/i(R) = oo, (2.2) 

which is clearly necessary for Z t to possess a distribution density. 
Let 

L(x) := x 2 l N<1 , if) L (0 := [ L^u)fi(du) = [ (Cuf^du), (2.3) 



'\€u\<l 

U(x):=x 2 Al, ip u (0-= I U(£u)n(du)= I ({£u) 2 A l) /i(du). (2.4) 



E 



Note that up to a constant multiplier the functions t/' 1 ' and ip u provide, respectively, a lower and 
an upper bound for the real part of the characteristic exponent 



Re -0(C) = / fl - cos{£u)) fi(du) 



(1 — cosl)< < -, \x\ < 1, 0<1 — cosx < 2, a; G 



Indeed, from the elementary inequalities 

1 — COS 2 1 

~x~ 2 - 2' 
we have 

(l-cosl)^ L (0<ReV(0<2^(0, (GR. (2-5) 
Our main assumption formulated below states that these bounds are in a sense exact. 

A. There exists some /3 > 1 such that 

^(0</^ L (0, ^K- (2-6) 

For £ > 0, we put 

Pt := (Re^) [ - 1] (1/t) = inf{^ > : ReV(0 > V0- (2-7) 



Note that under condition A the set in the right hand side of (2.7) is not empty (cf. Lemma 3.1 
below). 

Theorem 2.1. Suppose that condition A is satisfied. Then for every t > the distribution 
density p t of Z t exists, is infinitely differentiable, and for every to > there exist positive constants 
c = c(t ),d = d(to) such that (1.1) holds true. 

Remark. We will see in the proof below that Pt{ ) vanishes at ±oo, which means that the 
maximum in (1.1) is well defined. 
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2.2 Off-diagonal bounds 



Apart from the on-diagonal bounds of the form (1.1 ), one is typically interested in the off-diagonal 



bounds, which would control the behavior of pt(x) w.r.t. the space variable x. Below we propose 
a new type of such bounds, which we call compound kernel estimates. Let us discuss this topic in 
details. 

A typical form of an upper estimate one may think about, is the "bell-like" estimate 



Pt{x) < f t (x) 



(21 



with some kernel ft(x) which, when considered as a function of x for a fixed t > 0, is symmetric 
on K and decreasing on IR + . In addition, one would expect the kernel ft{x) to have the form 



ft(x) = a t f(a t x), 



(2.9) 



where the "shape function" f(x) is re-scaled by some "scale function" o t . Typical examples which 



give a strong motivation for the estimate (2.8) with ft(x) of the form (2.9), are the Brownian 
motion and, more generally, the a-stable process. In this case, (2.8) and (2.9) hold true with 
t and a t = t^ 1 ^, respectively. However, for more general Levy processes it may happen 



that any upper estimate of the form (2.8) with ft(x) as in (2.9), is very inexact. To see that, we 



recall the papers |Ish94j . |Pi97a| . [Pi97b], |Ish01] , where the pointwise asymptotic behavior of the 
transition probability densities of Levy driven Markov processes is studied. In our setting, the 
results from [Ish94j, |Pi97a]. |Pi97bj . |Ish01j imply that for some (not necessarily all) points x G 1R 
and for every to > there exist positive constants c = c(x,to), d = d(x,to) such that 

ct 7(x) <p t {x) < df( x \ te{0,t ]. 

Here the power 7 (re) may vary from point to point, which does not correspond well to the "bell- 
like" estimates ( |2.8 ), (2.9). In Example 4.2 below we extend this observation and give an example 
of the Levy process such that, if one constructs any estimate of the form (2.8) with the function 
ft{x) monotonous on M + w.r.t. x, then this function ft{x) is necessarily non-integrable on R for 
every t > 0. Since 

Pt(x) dx = 1, t > 0, 



we see that for this Levy process any estimate of the form (2.8) with a "bell-like function" ft(x) 



(not necessarily of the form (2.9)) is extremely inexact. Note that in Example 4.2, as well as in 
[Ish94], |Pi97aj . [Pi97bJ, |Ish01] . the Levy measure of the noise is a- stable like in the sense that 



ce 2 - a < 



u 2 fi(du) < de 



2-a 



ee (0,1]. 



(2.10) 



\U\<£ 



This means that even in the case, which is very close to the a-stable one, in general one can not 



expect to obtain a sufficiently exact upper estimate of the form (2.8) with a "bell-like function 



ft(x) without further restrictions on the Levy measure of the noise. 

To explain the proper structure of an upper estimate, we represent Z t as the sum 

Zt = Zf + Z t — a t , 



(2.11) 



where Z t and Z t are independent, Z t possessing for each t > the infinitely divisible distribution 
with the characteristic exponent 



M0=t I (l-e* u + t£u)fi(du), 

'\ptu\<\ 



(2.12) 
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Z t having the compound Poisson distribution with the intensity measure 



A t (du) = tp(du)l\ ptU \ >1 , (2.13) 

and 



a t = t I a + u p(du) J 

V Jl/pt<\u\<l J 



where p t is defined in (2.7). We show below that if condition A is satisfied, then Z t has a 



distribution density Pt(x) which satisfies (2.8), (2.9) with o t = p t . On the other hand 



p t (x) = (pt*P t *6-a t )(x) (2.14) 

with 

P t (dy) = e- A ^Y,- } K m (dy), (2.15) 



ml 

m=0 



where A* t m denotes the m-fold convolution power of the measure A t , and Ajf° is equal to the 8- 
measure at 0. This observation leads to an upper estimate for p t (x + at) in the form of the series 



of convolutions of the single kernel function ft(x) involved into (2.8) with the convolution powers 



of the measure A t . The above construction motivates the following general definition. 

Definition 2.1. Let a, £ : (0, oo) — > R, h : R — > R be some functions, and (Qt)t>o be a family of 
finite measures on the Borel a-algebra in R. We say that a real-valued function g defined on a set 
A C (0, oo) x R satisfies the upper compound kernel estimate with parameters (at, h, (t, Qt), if 

9t(x) <y2— } [ <Tth{x( t - y)Q* t m (dy), (t,x) e A. (2.16) 

Tfl\ /TED 

m=0 >/M 



If the analogue of (2.16) holds true with the sign > instead of <, then we say that the function g 



satisfies the lower compound kernel estimate with parameters (o~t,h,Q,Qt). 

In the above definition, the function h heuristically determines the shape of a "single principal 
kernel" involved in the estimate, which is re-scaled by the functions a t and ( t and then extended 
in a "compound Poisson" like way by the family of "transition measures" Q t , t > 0. 

Now we are ready to formulate the main result of this section. Denote 

Xt := min < x : Pt(x) = m&xp t (x'] 
[ v ' x'eM 

Below we will see that x t is well defined, and for every to > there exists L = L(t ) such that 

\x t \<L/pt t 6(0, t ]. 

Theorem 2.2. Let condition A hold true. Then for every t > there exist constants b±, 62, &3, W, 
such that the following holds true. 

I. The function 

p t {x + a t ), (t, x) G (0, t ] x R, 
satisfies the upper compound kernel estimate with parameters {pt, f upper, Pt, At) , where 

fu PP er(x) = b 1 e- b ^. (2.17) 



5 



77. The function 

Pt(x+a t -x t ), (t, x) G (0, t Q ] x R, 
satisfies the lower compound kernel estimate with parameters (pt, flower, Pt,-^-t), where 

flower(x) = ht\ x \<b 4 . (2.18) 

Under some additional restrictions on the Levy measure /x, an upper compound kernel estimate 
from Theorem 2.2 may lead to the "bell-like" estimate of the form (2.8). We formulate two 
statements of such a type in Theorem |2.5| below. However, as we have mentioned it in the 



discussion above, any statement of such a kind would necessarily be more sensitive w.r.t. to the 
structure of the Levy measure of the noise than that of Theorem 2.2[ which involves only our main 
condition A. This is the reason for us to focus on the compound kernel estimates. To finalize the 
discussion on that subject, let us mention that compound kernel estimates preserve some principal 
features which the "bell-like" estimates exhibit in the diffusion case. In particular, it is possible to 



extend compound kernel estimates, obtained in Theorem 22_ for Levy processes, to the Levy-type 
setting using an analogue of the parametrix method, see the forthcoming paper |KK12j . In the 
same paper we give some applications of the compound kernel estimates to the investigation of 
local times for Levy-type Markov processes. 



2.3 Extensions and corollaries 



In this section we formulate some extensions and corollaries of the above results. The first extension 
gives compound kernel estimates for the derivatives of the distribution density Pt{x). 

Theorem 2.3. Let condition A hold true. Then for every k > 1, to > 0, there exist constants 
61,62, such that the function 

gk 



dx 



: p t (x + a t ) 



, (t,x) G (0,to] x R, 



satisfies the upper compound kernel estimate with parameters (g t , f upper , g t ,A t ), where the func- 
tion f U pp er is defined in (2.11). 



The second extension gives more precise versions of Theorems |2.2| and |2.3| in the case where 
the process Z t is symmetric. 

Theorem 2.4. Let the process Z t be symmetric and condition A hold true. Then the statements 
of Theorem 2.2 and Theorem 2.3 hold with at and Xt replaced by zero, and f upP er defined by 

fu PP er(x)=b ie ~ b ^ 1+ W (2.19) 



instead of (2. 11). 

Finally, we give a corollary which provides the "bell-like" estimate for the distribution density 
Pt(x) and its derivatives. To do that, we recall the notion of a sub-exponential distribution, see 
|EGV79j and [Kl89j . 

Definition 2.2. A distribution function G on [0, 00) is called sub- exponential, if 
(i) for every y G R one has lim^^ = 1; 

fin ii m i-g* 2 (s) _ 2 

A distribution density g on [0, 00) is called sub- exponential, if it is positive on [x , 00) for some 
xo > 0, and 

(i) for every y G R one has linx^oo 9 ^ x ~^ = 1; 



(ii) lim, 
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Note that the tail of a sub-exponential distribution (respectively, a sub-exponential density) 
decays slower than any exponential function, i.e. 



lim V(x)e ax = oo, 



where V(x) = 1 — G{x) (resp., V(x) = g{x)), and a > is arbitrary. This follows from the 
representation 

V(x) = c(x)e~^ z o b ^ dy for x > zq, 



where zq = in the case V(x) = 1 — G(x), and zq = xq taken as in Definition 2.2 for V(x) = g(x). 



The functions c : [zq, oo) — > (0, oo) and b : [z , oo) — > R are such that lim^^oo c(x) = c > 0, 
lim^oo b(x) = 0, cf. |KI55| . 



Theorem 2.5. Let condition A hold true. 



I. Suppose that there exist a sub- exponential distribution function G and a constant C such 
that 



: \ptu\ > v}^j < C(l - G(v)), ve[l,oo), te(o,t }. 



(2.20) 



Then for every t > there exist positive constants bi,b2,Ci, such that 

Pt(x + a t ) <C lPt (f upper ( Pt x) + 1 - G(\xp t \)) , xeR, t e (0,t ], (2.21) 



where f U p per is defined by (2.17) 



77. Suppose that the Levy measure \i possesses the density m(u) 



fi(du) 
du 



, and there exist a 



sub- exponential distribution density g and positive constants C, to such that 



m t {u) := tp t m(p t u) < Cg(\u\), \u\ > 1, t e (0,t ]. 



(2.22) 



Then there exist positive constants b\, b<i, Ci, such that 



Pt(x) < G x pt (fu PP er(xpt) + g{\xp t \)) , x e R, t E (0,t ], 



(2.23) 



where f upper is defined by (2.17). 



Remark 2.6. For Theorem 2.5, the extensions analogous to those given in Theorem 2.3 and 
Theorem 2.4 are valid. Namely, for every k > 1, under the proper choice of the constants b\, b^, C^, 
we have 

<C k p$ +1 f(p t x), xER, te(0,t o ], ( 2 - 24 ) 



dx k 



Pt{x + a t/ 



with 



f up per(x) + 1 — G(\x\), when p satisfies (2.20), 
fu PP er(x) + g(\x\), when p satisfies (2.22). 



(2.25) 



For a symmetric process Z t we have a t = 0, and can take f upper in (2.25) defined by (2.19). 



7 



3 Proofs 



3.1 Auxiliaries 

We begin with a brief discussion of the properties of the auxiliary functions ip L and ip u . Everywhere 
below it is assumed that condition A holds true. 

Clearly, the functions L and U, involved in the definition of ip L and ip u , respectively, satisfy 

rxa 2 



U(x 2 ) - U( Xl ) 



X 



-L(x) dx, X\ < x 2 . 



Then, by the Fubini theorem, we have the following relation for ip L and ip u : 

<2 



hi 



(3.1) 



Combining this relation with condition A, we obtain the lower growth bound for the real part of 
the characteristic exponent ip. 

Lemma 3.1. Suppose that condition A holds true. Then there exists a constant c > such that 

Re^(0>c|£| 2//3 , |el>l- (3.2) 

Remark 3.2. By the standard properties of the Fourier transform, it follows immediately from 
(3.2) that for every t > the distribution density p t of Z t exists and p t G C£°(R). 

Proof of the lemma. Define 

9 u (v) := *p U (e v ), 9 L {v) := ^ L {e v ), v E R. (3.3) 



It follows from (3.1) that these functions are related as follows: 

fV2 



6 u (v 2 ) - B v {v l ) = 2 
Combined with condition A, this gives 



6 (v)dv, V\ < v 2 . 



e u (v 2 )-e u ( Vl )> 



6 U (v)dv, vi < v 2 . 



(3.4) 



(3.5) 



This means that for the Stieltjes measure on M + , generated by the non-decreasing function 9 , the 
absolutely continuous part in its Lebesgue decomposition is bounded from below by (2/ f3)8 u (u) du. 
Using this observation, one can easily show the inequality 



e -(2/0«20tf( Va ) > e-W^O u (v x ), Vl < v 2 . 
It follows from (2.2) and the definitions of the functions ip u and 6 U that 

e u (o) = ^(i) > o. 



(3.6) 



Therefore by (3.6) we get 



6 u {v) > Cl e (2/ ^, v > 0, 



with ci = 9 (0) > 0, which in turn implies 



Using condition A and the left inequality in (2.5), we get (3.2) with c = Ci(l — cosl)/3 1 . 



□ 



Together with the function p t defined by (2.7), we consider the functions 

p u t := inf{£ > : ij> u (£) > 1/t}, p\ := inf{£ > : ^ L (0 > l/t}. 



(3.7) 



In view of (3.2), (2.5) and condition A, the sets in the above definition are non-empty. Note that 

1 



iteV(ft) = V> u fa w ) = ^w) 



t 



(3.8) 



because the functions if) and i\F are continuous, and the function ip L is right continuous on 
In the sequel, we write f t x g t , t G A, if there exist positive constants ci, C2 such that 

ci/t <9t< c 2 ft, t G A. 
Lemma 3.3. For every c G (0, oo) one has 

p t ^ Pct ^p u t x ft L , te(0,l]. 
Proof. From the relation < ^ and condition A we have 

p U t <pf< p%. 



In addition, by (2.5) we have p^ t < pt < P(i- C osi)r Therefore, it is enough to prove that 



t G (0,1]. 



(3.9) 



Note that is monotone by construction. Hence one of the bounds in (3.9) is trivial (the lower 
one if c > 1, and the upper one if c < 1). The other bound in (3.9) follows from the inequality 



> 



2/(3 



6 > & > o, 



which holds true due to (3.6). 



□ 



3.2 Upper bounds 



Consider the decomposition (2.11). Note that the real part of the characteristic exponent ip t of 
the variable Z t satisfies 



ReVt(f) = tReij(0 -t (l - cos(£u)) n(du) 

J\p t u\>l V 1 



> *ReV(0 - 2t / /i(rfu) > tRe 
J| Pt «l>i 



1>(0-2til> u (jH 



Hence, by Lemma 3.1 for every t > the variable Z t has a distribution density G C, 



oo / 



dx k 



p t {X) 



2tt J r 



— j ?\ k P - ix t-'<l>t(0 



d£, k>0. 



(3.10) 

), and 
(3.11) 



By (2.14), the upper bounds in Theorem 2.1, Theorem 2.2, and Theorem 2.3 are provided by the 
following lemma, which gives respectively the upper bounds for the distribution density p t and its 
derivatives. 
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Lemma 3.4. For any k > 0, t > 0, there exist positive constants b\ and b 2 such that 

' Pt(x + a t ) < V" +1 e" 62PtN , XER, tG(0,t ]. 



dx k 



(3.12) 



Proof. We apply the complex analysis technique, similar to the one used in (KSTOaJ and |KK11] . 
Observe that the function ipt{0 has an analytic extension to the complex plane, and for given x 
and t the function 



H (t, x, z) = -ixz - ip t (z) = —i xz -t (1 - e luz + 

J\p t u\<l 



iuz)/j l (du) 



is well defined for z G C. We consider the integral in the right hand side of (3.11 ) as the integral 
of the function 

{-iz) k e H(t ^ z \ z G C 

along the real line in C, and change the integration contour as follows. 

Fix 7] G R such that \r]\ < pt, and for M > consider the rectangular contour T M = U*=i 
with 

rf = {y + iO, \y\<M}, r? = {y + i V , \y\<M}, = {±M + isrj, s G [0, 1]}, 

properly oriented. By the Cauchy theorem we have 

/ + / + / + / (-iz) k e H ^ x ' z) dz= [ (-iz) k e H{t ^ z) dz = 0. (3.13) 
/rf Jr!f Jrf Jvf\ Jr* 1 

Let us show that the integrals along j = 3,4, vanish as M — > oo. Consider the case j = 3, 
the case j = 4 is completely analogous. We have 



ReH(t,x,M + isr)) = sxr) — t / (l — usr] — e US71 )fj,(du 

J\ptu\<\ V J 



- 1 



\ptu\<\ 

e~ sriu {l - cos(Mu))n(du) 



(3.14) 



'\ptu\<l 

< sxr] — ip t (isr]) — e~ Re0t(M). 



Denote 



%(r]) :=if)t(ir]), C(t,x,r}) := sup (sxr]-%(sr))). 

se[o,i] v ' 



Then by (3.14) we have 



_ iz) k e H{t,x,z) dz 



1 3 



< (rf + M 2 ) fe/2 exp (c(t, x, 77) - e^Re^Af)) ■ 



By Lemma 3.1 and (3.10), the above inequality implies that the integral along T^ 1 vanishes as 
M — >■ 00. 



Passing to the limit in (3.13) as M — > 00, we get the representation 



(-iz) k e H{t > x ' z) dz = — j (-iy + v )k e vx-ix V -My+iv) dy _ 
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Thus, by (3.14) with y instead of M and s — 1, we have 



Or'-"' ' -2k 
Recall that \rj\ < pt, and denote 



(3.15) 



c := sup 



B 



1 — x — e 



\x\<l x 



Then, using the first inequality in (2.5) and the first identity in (3.8), we get 
-%(rj) <ct I {rjuf^du) < ctij L (p t ) < ( °—- j tRe^(p t ) = 



\ptu\<l 



1 — COS 1 



1 — COS 1 



Hence, we can take r\ = —p t sigax, and obtain from (3.15) 

d k 



dx k 



m 



(x) < Cl e- ptM / (p k t +\y\ k )e- e ~ lKe ^ y) di 
Jn 



with 



ifc-2 



Cl 



7T 



cxp 



c 



1 — COS 1 



Recall that Re^f satisfies (3.10) and note that, by (2.5) and condition A, the term 

t^ U (p t ), tG(0,l], 



involved in the right hand side of (3.10), is bounded. In addition, we have by (2.5) and condition A 
the inequality 

Rem > ^pV(o. 

Using these observations we finally get the estimate 



Qk 

dx k 



p t (x) 



<c 2 e-^(p k I (t 



1 — cos 1 

W 



+ h[t 



1 — cos 1 

W 



where 



I k (t,X) : = / \y\ k e-^ U ^dy, k>0, 



□ 



and C2 > is some suitable constant. The rest of the proof follows from Lemma 3.5 below. 
Lemma 3.5. For every k > 0, A > 0, there exists a positive constant c = c(k, A) such that 

I k (t,X)<cpl te (0,1]. 

Proof. Clearly, the function p t , t G (0, 1], is bounded from below by some positive constant. Since 
the function 

hVt 



4°(t,A):= / \y\ k e-^ U Mdy, te(0,l), 



i 

-l 

is bounded, it remains to prove 

I+(t,\)<cp k , lr(t,X)<cp k te(0,l], 
11 



I+(t,A):=/ y k e-^ U ^dy, /*(f,A):=/ (-y^^d 



where 



We prove the required relation for l£(t, A), the case of I k (t, A) is completely analogous. 
Making the change of variables y = e v , we get 



see the proof of Lemma 3.1 for the definition of 9 and other auxiliary functions. Take Vt = log 



Since 9 U is non-negative, we have 

POL 

< (pfr 1 / 

Jo 

The function 9 U is non-decreasing, hence 



e (k+l)(v-v t )-t\(e u (v)-9 u (v t ))-t\0 u (vt) dv 

e ( k+ i)( v - vt ya(6 u (v)-e u (vt)) dv 



rv t -i 

,(k+i)(v-v t )-t\(e u (v)-e v (vt)) dv < J e (k+i)(v-v t ) dv < 1 

~ — k + 1 



On the other hand, we have from the last identity in (3.8) and the trivial relation ip L < ip u that 
t9 u (vt) > 1. Then, using the monotonicity of 9 U , condition A, and applying twice (3.4), we get 



t[9 u (v) - 9 U (v t )} = 2t I 9 L {r)dr 



ft 



> 2t/?- 1 / 9 u {r)dr = 2tf3~ 1 9 u {v t )(v - v t ) + Atf3~ 1 



9 L (s) dsdr 



Vt >> Vt 



>2t(3- 1 U (v t )(v-v t )+4t(3- 2 
>2r 1 (v-v t ) + 2p- 2 (v-v t ) 2 . 



9 u (s) ds dr 



Therefore, 



,(k+l)(v-vt)-tX6 u (v)+tXe u (v t ) dv < / e (k+l)w-(2\/t3)w-(2\/l3 2 )wi dw 



< oo. 



Combining the above inequalities and using Lemma 3.3, we obtain the required statement. □ 



In the symmetric case we get essentially the analogue of Lemma 3.4 but with a slightly different 
upper bound. Note that in this case we have a t = Xt = 0. 

Lemma 3.6. Let ip be real valued. Then for any k > 0, t > 0, there exist bi, 62 > such that 

Qk 



dx k 



Vt 



(x) < 6 lp fc+V b2ptNln(ptN+1) , xgR, te(o,t ]. 



(3.16) 



Proof. The proof is essentially the same as in the non-symmetric case. First, we change the 
integration contour in the same way as in the proof of Lemma 3.4 Unlike the non-symmetric 
case, while doing this we do not pose any additional assumptions on rj £ R. The functions 
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H(t,x,z), ^t(v) an d C(t,x,rj) have the same meaning as before, we only need to keep in mind 
that now ipt{0 admits the representation 



?p t (0= / (l-cos(£iO)/i(diO. 

J\ Pt u\<l 

We have 

Re H(t, x, M + irj) = xi] + t / (cosh(r/u) — l)p(du) — t / cosh(7yw)[l — cos(yu))p(du) 

J\p t u\<l J\ptu\<l 

< xr] - ip t (iv) - tyt(M). 

(3.17) 



Thus, 



(-iz) k e H(t ' x ' z) dz 



< ( V 2 + M 2 ) fc / 2 e c( ^^-^ (M) , i = 3, 4. 



Therefore, applying the Cauchy theorem to the same contour as in the proof of Lemma |3.4| and 
letting M — > oo, we derive 



2tt 



(3.18) 



E. 



By Lemma 3.5 we get 



(\v\ + \y\) k e-^ y) dy< Cl (\rj\ k p t + p k t +1 ), 



]R 



where c\ > is some constant, depending on k. Since the Levy measure p is symmetric, we have 



-iptiif]) = t / [cosh(?ytt) — l]//(cfoi) = t / {7]u) 2 d{r]u)p{d 

J\ptu\<l J\ptu\<l 

(r)u) 2 p{du) = t{r)/p t ) 2 #(r)/ Pt )ip L {p t ) 



'\ptu\<i 

< t#(v/Pt) 



\p t u\<l 



cosh(r] /p t 



where $(x) = x 2 [cosh x — 1], and we used that i? is even and strictly increasing on (0, oo). 
Thus, by Lemma [3. 5[ we obtain 



dx k 



— p t (x) < Cl e-^ +cosh ^\\r,\ k p t + p k+1 ' 
The upper bound follows by minimizing the right-hand side in r\ over all rj G R. 



< /-„ n fc+l ^ — x»;+cosh(?7/ pt )+fc ln(?7/ pt ) 



< ^ n fc+l^-x»?+C3(fc)cosh(77/pt) 



(3.19) 



□ 



3.3 Lower bounds 

To obtain the lower bound, observe that the upper bound for pt{x) implies the existence of L > 
such that 



Pt(x) dx > 1/2. 



\ptx\<L 
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Then 



1/2 < / p t (x) dx <— max p t (x) 

J\ Pt x\<L Pt x ^ 



(3.20) 



which proves the lower estimate in Theorem 2.1 Using additionally that < bp\ with 



some b > by Lemma |3.4[ we obtain by Taylor expansion argument 

" d 



Pt(x) > Pt(x t ) 



~, dx 



Pt(y)dy 



> iiPt - c 2 pl 



This together with (2.14) implies the statement II in Theorem 2.2, because condition A provides 



that the family A t (R),t G (0, 1] is bounded. 



3.4 Proof of Theorem 12.5 

I. Put 



H(du) -=e- c J2 



C n G* n {du) 



n=l 



nl 



By (2.20) we can rewrite the upper estimate for pt(x) as 

/oo 
(fupper(xp t ~U) + fupper(xp t + V,)) H(du). (3.21) 



Since G is by our assumption sub-exponential, we have by Theorem 3 from |EGV79| that the 
distribution function H is sub-exponential as well, and 



C. 



x->oo l — H{x) 
which implies for x large enough 

c 2 (l - H{x)) < 1 - G{x) < c 3 (l - H(x)), 
where C3 are some positive constants. Then for any positive function / vanishing at 00 we have 

POO 

f(u)H(du) = I (1 - H(u))df(u) + /(l) 
(l-G(u))df(u)+f(l) 



f(u)G{du). 



Thus, we can write 



/oo 
(fupper(xp t -U) + fupper(xp t + «)) G{du). (3.22) 

Let us estimate the right-hand side for x > 0, the case x < is completely analogous. For x > 
and v > we have f upper (x + v) < f U pper{ x )-> which gives 



fupper(% V)dG{V) fupperip^) • 
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Integration by parts gives us for the other term in (3.22) 



oo 

(x - v)dG(v) = - I G(v)df upper (x - v) 



upper 

Jl 

/oo 
(1 - G(v))\f upper {x - V)\dv + fupp er (x - 1) 

x-1 

(1 - G(x - v))\f' upper (v)\dv + f upper {x - 1), 

-oo 

(here we understand the derivative in the sense of the derivative of an absolutely continuous 
function). Since G is sub-exponential and \f' upper \ is integrable, 

r*- 1 i _ Q( x - v \ r 

lim / ^7 \ \fLper( V )\ dv = / \fLper( V )\ dv - 



^ _ QhjA upper v / i \ J upper 

Thus, for x large enough we have 

oo 

(1 - G{v))\f upper {x - v)\dv < c 5 (l - G(x)). 



Finally, by the comment after Definition 2.2 we get (2.21). 
II. Denote 

^ m* n {u) 

n=l 



By (2.22) 



where 



is sub-exponential as well, and 



h t {u) < h(u), 



ni 

n=l 



li m = a (3.23) 

- 9{y) 



x—>oo 



Therefore, 

Pt{x) < Ptfupper{xpt) + C X pt I (fupper&Pt ~ V,) + fupper{xpt + «)) g(y)dy. (3.24) 



By the comment after Definition 2.2, the function g decays as x — )■ oo slower than any exponential 
function, implying 



{f U pper{xpt ~y) + f upper (xp t + y)) g{y)dy <c 2 g{\xp t \). 



Thus, we can estimate the right-hand side of (3.24) as 



C 3 Pt ( J upper (xpt) +g(\xp t \)) 



which finally gives (2.23). □ 
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4 Examples 



Example 4.1. Let Z t be a symmetric a-stable process, a G (0,2). In this case t/>(£) = |£| a , 
p(du) = C a \u\~ a ~ 1 du, and one can easily verify that condition A is satisfied. 

In this case p t = Because the Levy measure possesses the density m(y) = C a \y\~ l ~ a , 

we are in the situation of Theorem 2.5 II with g(u) = (l/a:)|tt| _1_Q: l u >i, which is obviously a 
sub-exponential density. Applying this theorem, we arrive at the well-known upper estimate for 
the symmetric a-stable distribution density: 



Pt(x) < Cl r l ' a f e -M-v« N in( t -v« N+1) + (t-i/«\ x fj a iy j 



< c 2 t 



-1/a 



A 



In addition, it is straightforward to verify that the lower compound kernel estimate from Theorem 



2.2 now provides the similar lower bound 



Pt(x) > c 3 t 



-l/a 



A 



l+o 



with some positive constant C3. In other words, using our main Theorem |2.2 we can re-establish 
the well-known two-sided estimate for the symmetric a-stable distribution density: 



Pt 



(x) x r 1/a f(r 1/a x), (t, x) e (0, 00) x 



with 



f(x) = 1 A \x\- a -\ 

Example 4.2. (a) Consider a Levy process with the discrete Levy measure 



(4.1) 
(4.2) 



00 

H{dy)= 2 m (s 2 - n v(dy) + 5_ 2 -n V (dy)y 



where v > 0, < 7 < 2v. Straightforward calculation show that in this case (2.10) holds true 
with a = j/v. In particular, condition A is satisfied and p t x t~ l / a . This means that there exists 
q > such that for every t G (0, 1] the inequality 



n < n (t) :- 



1, 1 

- log 2 - - q 

.7 t 



implies 
Consequently, 



2- nv Pt > 1. 

At(dy)>t 2 n i(8 2 -n V (dy) + 5_ 2 -nv(dy) 



n<no(t) 



and, taking into account that in this case at = Xt = 0, one has by Theorem 2!1_ II that for some 
positive constant c 

Pt(2- nv ) > ctp t 2 n \ n < n (t). 



Then every function ft(x), monotonous with respect to x on R + and satisfying (2.8), should satisfy 

ft{x) > ctp t 2 n \ x G (2-™ 2-( ft - 1 )«], n < n (t). 
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It is easy to verify that when 7 < v (which is equivalent to a < 1) one has 



n<no(t) 



->• OO, t ->■ + . 



This gives J R ft(x)dx — > 00 as t — > 0+, which shows that in this case any estimate of the form 
(2.8) with a "bell-like function" ft{x) would be extremely inexact. 



(b) Consider now the case 1 < a < 2. Note that 

tJ{u : \ Pt u\ >v}\ =t T<n ^ Ct\A~ 1,v = Ct\x\~ a , x>2 v , 



n x <n<0 



where n x := — log 2 Therefore, condition (2.20) of Theorem 2.5 I holds true with 1 — G( 



[x 



\x\ a . By this theorem and Remark 2.6 , we have the following estimate for the respective transition 
probability density: 



with 



Pt{x) < c^l-^e-^- 1 ^^- 11 ^ 1 ) + (t- l ' a 

f(x) = 1 A \x\~ a . 



x 



-l/at 



X 



(4.3) 



(4.4) 



This upper bound looks very similar the one from (4.1 ), with the notable difference in the functions 



/ in (4.2) and (4.4). This difference appears because in the previous example we have used the 



version of Theorem |2.5| based on sub-exponential densities, while in the current example we had 
to use another version of this theorem which deals with sub-exponential distribution functions. 



Note that (4.3) is precise in the sense that 

Pt{x) > c 2 r 1/a f{t- 1/Q x), t>0, \x\ G {2- nv ,n > 1} 



(4.5) 



with / given by (4.4 ). The latter inequality can be verified easily using the lower bound in Theorem 



2.2 On the other hand, one can hardly expect that (4.5) can be extended to hold true for all 



iGK, which would be a complete analogue of the lower bound in (4.1) 



The third example shows that condition A is comparatively mild. In particular, it allows some 
kind of "slow oscillations" for the characteristic exponent. 



Example 4.3. Consider the function a : [0, 00) — y [a_,a+] C (0,2) such that 

va'{v) — > as v — > 00. 

Put 



e-(v) :-- 



e a{w)w dw, 6 + {v):=2 9~(w)dw. 

n JO 



Using the FHospital rule we get by (4.6) 



e a(v)v 

9~(v) ~ — -— as v — > 00. 
a(v ) 



(4.6) 
(4.7) 

(4.8) 



Then the function 0(£) := 9 (0 V ln£) satisfies 



a) 4,( < (which follows from §Ut\ and (pf); 



b) ^ as f -> 00. 
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Thus, there exists a Levy measure fi(du) = m(u)du such that 



0(0 



(C,u) 2 m(u)du, 



l€u|<l 



or, in other words, </>(£) = ip L (!;), where ip L (£,) is the lower function, corresponding to /i. Then 
the corresponding upper function is ip u (^) = 6+(0 V lnf). Similarly to (Ob we get 



which implies together with (4.8) 



9 + {v) 



V(0 



2 e a(v)v 

a 2 (v) ' 



v — > oo, 



^(0 a(ln0 
Thus, ^ satisfies condition A, and from (4.8) and (4.9) we get 

^(£)x^(0x^(0x|r (ln|fl) , |£|->oo. 
In particular, for any a G [a_, aq_] there exists a sequence {£ Q ,n, w > 1}, such that 

£ a , n ->■ OO, V(^a,n) X (^a,n) Q , 1400. 



(4.9) 



(4.10) 



Thus, by Theorem 2.1 for any a G there exists a sequence t a . n — > such that 

P^(0)x(t Q , n )" 1/a , n^oo. 

In this example, heuristically, the asymptotic behaviour of the characteristic exponent ip may 
coincide with the one of the characteristic exponents of a-stable processes on various subsets of 
R + having oo as their limit point. Here the index a depends on the subset, and may vary in the 
range [a_,a+]. Respectively, the asymptotic behaviour of the transition probability density as 
t — > 0+ depends heavily on the particular subset of R + this behavior is considered on, and may 
be the same as for the a-stable processes with any a from the segment [a_, a + ]. 



Acknowledgement. The authors thank N.Jacob and R. Schilling for their helpful comments 
and remarks. 
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